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In order to study and control the attitude of a spacecraft, it is necessary to understand 
the natural motion of a body in orbit. Assuming a spacecraft to be a rigid body, dynamics 
describes the complete motion of the vehicle by the translational and rotational motion of the 
body. The Simulink Attitude Analysis Model applies the equations of rigid body motion to 
the study of a spacecraft’s attitude in orbit. Using a TCP/IP connection, Matlab reads the 
values of the Remote Manipulator System (RMS) hand controllers and passes them to 
Simulink as specified torque and impulse profiles. Simulink then uses the governing 
kinematic and dynamic equations of a rigid body in low earth orbit (LEO) to plot the attitude 
response of a spacecraft for five seconds given known applied torques and impulses, and 
constant principal moments of inertia. 


I. Introduction 

Considering a rigid body, defined as a system of particles where the relative distances between particles is fixed, 
French Mathematician Michel Chasles states that the body’s translational motion along a line and its rotational 
motion about that line completely describes its motion. Therefore, the orientation of a rigid body is described by its 
position, velocity, angular velocity, and angular position. The angular velocity and position of a body is further 
described by Leonhard Euler’s theory of rotation, which states that the rotation of a rigid body may be described by 
the three Euler angles. Applying Chasles’ Theorem and Euler’s Rotation Theorem to the motion of a spacecraft in 
Low Earth Orbit (LEO), the Simulink Attitude Analysis Model reads the values from the RMS hand controllers 
using a TCP/IP connection and determines the vehicle’s final position and velocity given known applied impulses 
and models the attitude response of the vehicle given known applied torques. 


1 Fall Intern, EV82 Stage Analysis, Marshall Space Flight Center, University of Alabama in Huntsville 
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II. Rotational Motion 


A. Using Matlab 

According to Euler’s Rotation 
Theorem, any rotation can by described 
by three Euler rotation angles. There are 
several sequences of Euler angles that 
depend upon the axes of rotation. The 
Simulink Attitude Response Model uses 
the most common Euler angle sequence 
known as the x-convention sequence or 
more commonly as the Euler 3-1-3 
sequence. As seen in Figure 1, the first 
rotation is by the spin angle (0) about the 
z-axis, the second rotation is performed 





Figure 1: Euler 3-1-3 Sequence 


function aodot =rotat ionaldynaniicE [t t x] 
1x^4 00; %(kg*iQ^2} 

Iy=50Q f %tkg*ra~2} 

% (^*111^2) 

T= [1000 ;0;- 1000 ] j % 
xdot (1, 1 )=jc( 2) /Ix+TU 

xdot (2 r l)=x;U) *x(U fly+ T[2 

xdat ( 3 F 1) =x < 1) *x C 2 ) * Elx-Iy) /Iz+T[3 
xdot (4,1)= (a in (x(6) ) * x ( 1 ) +cos [x ( £ ) 
xdot (5 r l)=oosCx{6) ) # xtl) -sin [x<6 ] ) 
xdot (E, 1) =x (2) - {sintx EG) ) *cos Ex (5 J 


/lx; 

/Iy; 

/Iz; 

*x ( 2) } / &in (x [ 5 ) 3 ; 
x[2} ; 

*x(13 +cos(x [£ J ) *coa (x ( 5] } *x (2) ) /sin (x [5) } 


by the nutation angle (0) about the former x-axis, and the third and final rotation is by the precession angle (\p) about 
the former z-axis. 

The function 
developed in Matlab 
to simulate the 
attitude response of 
a spacecraft to 
known torques is 
seen to the right in 
Figure 2. The 
program 
rotationaldynamic s . 

m provides the Figure 2: Rotationaldynamics.m 

governing differential equations of motion. The first three equations, xdot(l,l), xdot(2,l), and xdot(3,l) calculate 
the angular velocity of a spacecraft given a known torque profile. The last three equations, xdot(4,l), xdot(5,l), and 
xdot(6,l) calculate the response of the three Euler angles, 0, 0, \|/ to the same known torque profile. Then using a 
variable step Runge-Kutta method, Matlab ’s built in 
ode45 solver numerically integrates the six 
governing equations. The results are seen to the 
right in Figure 3. 

B. Using Simulink 

Developed by Mathworks as an add-on product 
to Matlab, Simulink is a graphical block- 
diagramming tool used to simulate dynamical 
systems. It provides a comprehensive library of 
block diagrams that the user drags into a workspace 
to build graphical diagrams as opposed to Matlab ’s 
numerical computing and programming 
environment. In order to integrate the rotational 
equations developed in Matlab to a graphical model 
in Simulink, the Attitude Response Model utilizes 
the S-Function block. S-Functions are useful 

because of their ability to incorporate existing C code into a Simulink diagram. Similar to rotationaldynamics.m that 
was built in Matlab, newrotationaldyn_sfun.m, seen in Figure 4, describes the system as set of mathematical 
equations. Using a set of callback methods, the S-Function block performs tasks at various stages of the simulation. 
The first step required is the initialization task, which sets the number and dimensions of the inputs and outputs, 
establishes the sample time, and allocates necessary storage. Then the S-Function calculates the next sample value 
and the output of each time step. Finally, the S-Function integrates the functions in the same way that Matlab’s 



Figure 3: Matlab Results 
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ODE solvers do. After writing an S-Function file, the inputs are defined in the Simulink environment and necessary 
parameters are passed to the file. Using the scope block, the results are displayed as x-y graphs. 


function [eya p xQ p Etr , Ls] - newxotat ionaldyn af (in £t T x p u e flag, xO) 
switch flag 

case 0 % initialization; 

sizes = slmslzes; 

slzos .KumCont State* = S; 

siz&a . UunsDiscStat-os - Q ; 

s 1 zes . SlurciQutput a - £ ; 

sizes j U mnInputB 

si z,es . Dir Feedthrough - Oj 

s i zes . NumSampleT imes = 1 j 

sys = siTT.si zes [s izes 3 ; 

scr = cl i 

te W [O 01 ; 

case 1 % derivatives; 

% It is convenient to use camion notation for states 

wx-x ( 11 ; 

wy= x ( 2 ) i 

wz=t(3} j 

pSl-=X(4) j 

Lheta=xf 5 > f 

pm=x Ct] ; 

% also, the inputs aro 

Inertia - u(l-3J ; 
roll^-u (4 ) f 
pitch^u{5) j 

yaw=u(6 ) i 


% state derivatives 

xdotl- (wy*wz* (Inertia (2,13 -Inertia (l p 1) J /Inertia ( 1, 1) *roll/lMrtia (1 „ 1) 3 *10O/pi; 
xdot2- (wx*wz* £ Inertia (2 , 13 - Inertia ( L p 1) ) /Inertia £2 , 1 ) + pitch/ Inertia (2, 1 J ) * 190 /pi } 
xdotl- (wx*wy* £ Inertia (I „ 13 Inertia (2 p 13 3 /Inert ia ( 2 , D * yaw/ Inertia ( 3 1 J 3 * IB 0/pi ; 
xl= { ( ain (phi) 'ra+coa (phi ) *wy) /sin (theta) ) * 100/pi ; 
x2= (cos {phi 3 *wx-sin (phi 3 *wy3 *l&Q/p!j- 

x3-(wz- (sin {phi 3 *cos ( theta] [phi 3 *cos (theta) *wy] /sincthata) }* 180/pi E 
sys = [Xdotl XOOt 3;Xl ;x2;x3] ; 

case 2 It outputs; 

eys - [x(l)jx(2) ;x{3J ,x(4> ; x ( 5) ; x (£ > j , 

case | 2, 4 r 9} 
ays = [] i 

otherwise 

error ( [ " unhandled flag = 1 , nuiri2str {flag) ] ) ; 
ond 


Figure 4: S-Function File 


II. Translational Motion 


A. Using Matlab 

Disregarding the rotational motion of a rigid body and focusing solely on the flight path of the vehicle allows for 
analysis to be confined to the translational motion of the spacecraft in orbit. When studying translational motion it 
is important to reduce the vehicle dynamics to the motion of a specific point on the body and study the entire body 
as a point mass system. The position, velocity, and acceleration of the particle represent the three degrees of 
freedom of translational motion. 

The translational function written in Matlab explores the governing equations that rely on the rotation matrix and 
four orbital parameters. The first step is to derive the angular momentum vector given the initial position and 
velocity of the body in motion. After determining the magnitude of the velocity impulse, which depends on the 
norm of the initial velocity and the orbit inclination, the velocity change can then be determined. The final velocity 
vector is then calculated by vectorially adding the velocity impulse to the initial velocity. The translation.m file and 
the rotation.m file are seen below in Figures 5 and 6. 
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clear, elej 

i:oM'. = r.l fgf sank 

riM-iMiDiMb mum :vb- 

v$»ifi-3rlM3T|9 J f liulttat wdocity ica-s; 
to 10 ; W£-g, orbit li wlt Tiatlffl 
gsrlZD; VdM, :»gti iscu^siar. ot isracdlag nrf a 
¥-Ib,‘ kdeg, j :i parlapais 

botixlODi ldflg 4 JDgia it tizich vnlccity iapil® li ippliirt 
dtr-pi/l&D; 

Ci - notation ! i 'd: c, cm' Itr, v'dcr] ; 

rav-Ci'rilj 

vafl-Ci^vd; 

hi=CTC«s ■' ra’3^ VaO^ - iaiyTiilar uoraeTitTiD dC Initial orbit Ikn/s « 
virs= 12222 4444 :8^] ■ 

Mipdr. =EKHL';tin) r Inaflnltude ul TSiccity itjkjIee e! 

frnfi?# 1 Iddb Ibeta'dtrl *TaC 7 , ®nUYaC , )+Bi:i(be':,a*dtr) 'lU/DQiafbil ) f lYtlcCity ckangE 
Yf^vaD-dV Uiaal relGCity !Wfil 


?f T 

a.msj 

■ILOSlf 

1.HS9 


fiiTicctwi C^matiwHl.dn.nj 

IDefines the rotation uatrlj of 3-L-j Biler angles [radians) 

n^er-s l=ort>ltal inclination.. om=i:lg2it ascension of ascending nods, and 

lv=argimieiit of periapsls 

LI-tos (an) "c&s I v> -sin l¥l " cos ( 1 ] ; 

[i2u -cos (cm) *s in \v\ sin ■! cm) * nos (v) *ooa [i ) j 
L3=sin(an) p sinrih 

>11 - sin (cm) *coa [*!■ i cna I nn.] * Ein (v) *coa [ i ) ; 

M2 = - sm (on) *s 1:1 \*\ #cos ( c*o * cos (v) ‘cos (ii ; 

>13- -cos (cm) *ain[il f 
B l=sin(¥)**iiiU1ij 
h 

H3 = Wfi(l); 

C=[L1 U L3;Ki N2 KjBl H2 Ml]; 


IqpuE argument “cm" is undefined. 

Error in ■*> rotation at 5 
tl=cos ftMj -fin rvP *ccs fij; 


Mfchrfiriiti J&ilim 'JO 


Figure 5: Translation.m Figure 6: Rotation.m 


B. Using Simulink 

As seen below in Figure 7, the translational motion modeled in Simulink uses constant blocks to generate inputs 
and mathematical blocks to modify the inputs. Lines are used to pass the values into the embedded Matlab function 
block as inputs. The embedded Matlab function allows the user to write a Matlab function that is used in a Simulink 
model. 





Figure 7: Translational Motion in Simulinnk 


III. Connecting Simulink to RMS Hand Controllers 

TCP/IP (Transmission Control Protocol/Internet Protocol) is the basic communication language of the Internet. 
TCP/IP allows for services such as file transfer across a very large system of servers. TCP converts files into data 
packets that are transmitted over the network connection to the destination computer. The IP simultaneously 
configures the destination address to verify the data is sent to the right place. In order to open a TCP/IP connection 
and read the values of the RMS hand controllers, the program seen in Figure 8 was written. It requests a service 
from a server in the network that reads the request, processes the data, and sends it back to the computer in a format 
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that Matlab can understand. The values are then passed to the Matlab workspace where Simulink reads them as 
torque and impulse inputs. 


kjbalizjc Lhf variables byte = uifllSCO), 

t initialize the file ■descriptor 
fd = ccpiprnns.-msfc.nasa.gov 1 , GSQQt ; 

I set the buffer size 

set (fd, 1 InputEufferSize 1 , 37); 

\ open, the connection 
fopen tfd) ; 

I tell the server to send the first record 
\ f print f (fd, 'u ']; 


Error using -■> icinterface'.fopeiz at S3 
[ftifcjiowsi flemolteii’osC : nss . msf c. na.sa . gcv . 

Error in ==> rinsnew at 13 
fopan fftfj 


% read some records 
for I = 1:100 
dispM) i 

\ tell the server to send a record 
fprintf [fd. "u 1 J : 

V read the record 

Idata, bytesl = f scant ff d, 1 Its r , 15] «■ 
dispidataJ e 

s = data ( \ ) ; 

nl - mod( ninte (data f 2 \ 1 ,641 r 
n2 - mod* uintB (data O > ) ,641 r 
nl - modi! ulntB (data [i|)j ,641 r 
nl - mod( uintB (data [5} ) ,641 e 

x =. inti 6 (nil* 40 S£ * iatl6 [n2] *256 * intl£ (nl] *16 * intl6 [n4] ; 
if [n i» '■'} 
x = x * -If 

etidj 

s = data[S); 

nl = modi uintB (data [7]) r 643 1 
n 2 - mod) uintB (data [0]) r 64] j 
n3 - noli uintB (data [9]) r 643 j 
:A = (rod; uintB (data [10) ) Mh 

y - IntlS (nl|MD9C + M16[n2m + MIS I nl) ‘IE + intlfi[n4) 

if (s == '-'I 

y = y * -lj 

end; 

s = data ( 111 ? 

nl = modi uintB (data [12)) r 64) ; 
n2 = modi uintB(data[13)l r 64) ; 
n3 = mod) uintB (data [14)) r 64) ; 
n4 = mod) uintB (data [15)) r 64); 

2 = IntlS (nil *4096 + lntl6 [n2J*256 + MlEtnai'lt + intlfi [n4] 
if (s ,= 
i = z 1 -If 
end? 

s = data 1 16); 

nl . modi uintB (data [17)) r 64); 


Figure 8: TCPIP.m 
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